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(57) Abstract: The invention is directed to a method for diagnosing the state of a system. The system may be mechanical, chemical, 
electrical, medical, industrial, business operations, manufacturing related, and/or processing related, among other. The method may 
measure a signal from the system. Further, the method may compare the signal to an expected signal. The method may calculate a 
signal strength and/or a noise. The signal strength and noise may be functions of a frequency. Further, the signal strength and noise 
may be used to determine a channel capacity and/or a rate of information. A comparison ol tlic rale of information and the channel 
capacity may yield information associated with the state of the system. The information may be used in diagnosing the state of the 
system. Further, the expected signal may be derived from a model. The model may be tuned to the measured signal. The model may 
have parameters that are associated with features and/or faults of the system. These parameters may be used in diagnosing the state 
of the system. Further, selectively repeated diagnosis over time may yield a prognosis of the system. 
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MODEL-BASED MACHINE DIAGNOSTICS AND PROGNOSTICS USING THEORY OF 
NOISE AND COMMUNICATIONS 

5 BACKGROUND OF THE INVENTION 

Related Applications 

This application claims priority of U.S. patent Application, Serial No.60/235 ,25 1 , filed 
September 25, 2001 entitled: "MODEL-BASED MACHINE DIAGNOSTICS AND PROGNOSTICS 
USING THEORY OF NOISE AND COMMUNICATIONS", and is incorporated herein by reference in 
10 its entirety. 

1. Field of the Invention 

The present invention generally relates to a method of diagnosing systems. In particular, the 
present invention relates to a model-based method for diagnosing the operational health of a system, 
15 and for forecasting the future operational health of the system 

2. Description of Prior Art 

Diagnostics and prognostics are used in many fields. These fields may include mechanical, 
chemical, electrical, medical, manufacturing, processing, and business operations, among others. Each 
of these fields has problems and difficulties relating to deterrnming the source of a problem, identifying 
20 the severity of the problem, and pretfcting the behavior of a system in relation to the problem 

For example, reliability and maintenance of complex equipment is critical to productivity and 
product quality. The purchase price of many typical equipments may account for half the equipments' 
cost Maintenance and support during the "life-cycle" may consume an amount roughly equal to the 
book value of the asset 

25 For example, at a typical chip plant billions of dollars are invested in equipment; many 

traditional manufacturing plants invest hundreds of millions of dollars in mamifacturing machinery. 
Most rnaintenance is rigidly scheduled by "time in service", not condition Machine productivity is lost 
during rnaintenance downtime. This and unscheduled downtime due to failures represent a very large 
part of a machine user's cost of operation. 

30 In these embodiments, machines are complex systems of components: gears, shafts, bearings, 

motors, lead screws, sensors, electronics, microprocessors, etc. integrated into a working whole. 
Machines in this context can also be biological, chemical, or hydraulic, among others. Defective or 
degraded components, alone or interacting, can render a machine dysfunctional. The machine may fail 
catastrophically and not complete its task, or it may lose tolerance, resulting in defective parts. 

35 Although methods and models exist for many individual component failures, errors from slightly 
degraded components can "stack", yielding overall system failure even when these models predict 
health of all individual components. 
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Many typical designers of machinery and engineers that maintain machinery, essentially know 

what the potential system faults are, and at what locations in the machine these faults will occur. 
However, unexpected breakdowns will still occur. Many typical designers and engineers do not know 
and/or can only poorly predict when &ults will occur. Further, periodically, healthy machinery must be 
5 taken out of service for maintenance. Thus adding an unnecessary cost Also, it may be very difficult 
and sometimes impossible to observe many of the conditions internal to the machine that lead to failure. 

One can easily imagine inetaphoric extensions of these problems to other fields such as 
chemical, electrical, medical, rnanuiacturing, processing, and business operations, among others. As 
such, many typical systems suffer from deficiencies in providing accurate diagnostics and prognostics. 
10 Many other problems and disadvantages of the prior art will become apparent to one skilled in the an 
after comparing such prior an with the present invention as described herein. 
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SUMMARY OF THE INVENTION 

Aspects of the invention are found in a method to diagnose operational health of a system, and 
to forecast future health. For example, the method may permit intelligent scheduling of mamtenance 
downtime in a mechanical or chemical system Further, the method may be used, for example, to avoid 
5 functional and catastrophic failures. 

Further aspect of the invention may be found in the method assembling models of the machine 
system, including system components and known system faults. Faults may be treated as "noise". In 
addition, parameters in the model may be "tuned" from signals from the real system, causing the model 
to mimic the real system in its present condition. Diagnosis may then be performed by observing the 
10 model 

In another aspect of the invention, the method may treat the system as a communications 
channel, estimate signal and noise levels, and diagnose health of the system with a tuned model by 
assessing how much information per unit time the system in its present condition can convey over its 
comrnunications channel. The method may compare mis maximum amount to the amount required by 

15 the system to execute a certain task. The method may assess if a system, in a given state, can perform a 
certain function, with a specified performance, within a given tolerance. 

Other aspects of the invention may be found in a method based on fundamental principles of 
physics and information theory. Further aspects may be found in the method assessing functional 
condition or state operable to perform a specified task, in addition to potential for catastrophic failure. 

20 Additionally, the method may operate on a tuned model, avoiding interpretation of complicated signals. 
Furthermore, the method may allow predictive scenarios for a system's possible future health and 
functional condition, given certain observed trends. 

In another aspect, for a different system or a new system design, only the model may be altered 
and not the basic diagnostic algorithm The method may also permit the incorporation of knowledge of 

25 faults, and the intent of the designers of systems, into the diagnostics routines. 

Another aspect of the invention may be found in assembling detailed dynamic systems models 
of the system in question- The models may posses a one to one correspondence between at least a 
portion of the physical components or elements in the real physical system and elements in the dynamic 
systems model. In some embodiments, all of the physical components are modeled in a one to one 

30 correspondence. 

In one embodiment the model may include all possible faults and potential failures in the 
system models. These effects may be tabulated as "noise" in the system Noise in a signal is the 
difference between the actual signal and the expected signal. In the model, noise may be induced by 
changes in parameters of dynamic system elements, which then alters any signals passing through a 
35 system Alternatively, if a certain fault cannot be described by these means, then sources of noise may 
be inserted into the system model, at locations in the model that are consistent with the locations of the 
faults in the real system. The intensity of these noise sources may be adjusted to make the model behave 
like the real system 
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In an exemplary embodiment, the method may former include, placing sensors on a machine, 

to monitor the machine; exciting the machine; and observing the machine's response via the sensor 
outputs. The collected data may be used to tune the model's parameters, so the model niimics the real 
system. After data has been collected on the actual machine, the system model may be excited with the 
same excitation as the actual machine. The outputs of the model may be compared to the corresponding 
outputs of the real machine. If the model's outputs differ from the real machine's outputs, the values of 
model parameters may be adjusted or changed, including the intensity of the noise sources, until the 
model's outputs approximate the actual system's outputs. 

In one embodiment, the channel capacity, C, of the system may be calculated The channel 
capacity may be the maximum amount of mformation per unit time that can be measured from 
successfully conveyed through the machine. The channel capacity may depend on the design and 
construction of the system, and the present condition of the system, which results from aspects. These 
aspects may include manufecruring, aging and damage, among others. For example, faults may be 
encoded as "noise" in the model. Analytically, the channel capacity may depend on the strength of the 
noise levels in the system, relative to the strength of me excitation system response signal. 

In the exemplary emrxxtiment, for a desired task to be performed by the machine, the rate of 
informaiion associated with the task may be calculated The rateof mformation may depend on the 
desired speed at which the machine does the task, the desired loads imposed on the rriachine, the 
complexity of the task, and the desired accuracy at which the machine should do the task. Further, the 
20 rate of information may be measured 

Another aspect of the invention may be found in comparison of the rate of information to the 
channel capacity. This comparison may be used to evaluate the operabihty of the system If the rate of 
ir*ormarion is less than or equal to the channel capacity, the system may perform the desired task within 
the desired precision. If the rate of information is greater than the channel capacity, the system may 
25 functionally fail. 

Another aspect of the invention may be found in the lonnulation of extremely ucuiiled models 
of the system to describe a system's condition. In one exemplary emtodiment of a system, the model 
includes bond graph based models of a motor, a gear box, and other mechanical transmission 
components. These extremely detailed models (a) exhibit a one to one correspondence between 

30 elements in the model and components in the real system; (b) incorporate many typical effects of the 
device into the model, including defects; (c) include in the models, via finite element concepts instilled 
into bond graphs, the dynamically distributed nature of components in the real system, and (d) use noise 
sources to account for defects and degradation of components. Simulation of the motor and gear box 
models may generate the complex spectra measured during operation of these devices. These models 

3 5 may rnirnic real system behavior and may be used to store information regarding the health condition of 
the machine. 
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In a further aspect, the models tabulate the effects of system faults (system maladies) as 

"noise" in the machine. Noise may be the difference between the actual signal received, and the 

expected signal that should be received As a machine degrades or ages, the difference between actual 

and expected signals may become larger. Thus noise levels may increase. These noise methods permit 
5 incorporation of faults into the models that heretofore could not be described analytically. The herein 

described methods have imported this body of knowledge to mechanical, hydraulic, other physical 

systems, and others, to name a few. 

In an additionally aspect, the method maybe used to predict the future conditions of systems, 

for scheduling maintenance and avoiding functional and catastrophic failures of the systems. The 
10 method may forecast if a complex system is capable of perfonning a given task, at a given speed and 

load, within a specified tolerance. 

The model system may be implemented on a computer system Hardware and software 

components may in combination allow the execution of computer programs associated with the method. 

The computer programs may be implemented in software, hardware, or a combination of software and 
15 hardware. 

Further modifications and alternative emfrxxnments of various aspects of the invention will be 
apparent to those skilled in the art in view of mis description. Accordingly, this description is to be 
construed as illustrative only and is for the purpose of teaching those skilled 

As such, a method for diagnosing and prognosticating the state of a system is described- Other 
20 aspects, advantages and novel features of the present invention will become apparent from the detailed 
description of the invention when considered in conjunction with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic block diagram depicting the Shannon-Weaver Model for use according 
to the invention. 

Figure 2 is a schematic block diagram depicting the informauon path according to the invention. 
5 Figure 3 is a schematic block diagram depicting a series of information paths according to the 

invention 

Figure 4 is a block schematic diagram depicting a computation system for imple menting the 
method, according to the invention. 

Figure 5 is a schematic block diagram depicting a network system for implementation of the 
1 0 msihod, according to the invention. 

Figure 6 is a block flow diagram depicting an exemplary method according to the invention. 

Figure 7 is a cross sectional view of squirrel cage induction motor. 

Figure 8 depicts Ghosh and Bhadra's [5] bond graph of a squirrel cage induction motor. 

Figure 9 depicts the stator resistances in Figure 8 redistributed to each of the stator coils. 
1 5 Figure 10'depicts a simplified representation of the signal and modulated GY element 

Figure 1 1 depicts a squirrel cage rotor with five bars. 

Figure 12 depicts a transformation of a and J3 phase currents into rotor bar currents. 
Figure 13 depicts the bond graph structure includin g stator and rotor bar action. 
Figure 14 depicts the bond graph equivalence used in modeling. 

,i» — i 

20 Figure 1 5 depicts the bond graph representing stator and rotor bar action in the magnetic circuit 

Figure 1 6 depicts angular velocity of rotor axis and stator currents in stator winding. 

Figure 17 depicts angular velocity of rotor axis and stator currents in stator windings, at startup. 

Figure 1 8 depicts angular velocity of rotor axis and 5-currents in each rotor bar, at startup. 

Figure 19 depicts angular velocity of rotor axis and 5-currents in each rotor bar, at startup. 

25 Figure 20 depicts angular velocity of rotor axis and 5-currents in each rotor bar, from startup to 

steady state. 

Figure 2 1 depicts stator currents and rotor velocity of a machine with a broken rotor bar. 
Figure 22 depicts stator current of 2nd phase and rotor velocity of a healthy machine at steady 

state. 
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Figure 23 depicts staler current of 2nd phase and rotor velocity of a machine with a broken rotor 
bar at steady state. 

Figure 24 deplete the angular velocity of rotor axis and 5 currents in each rotor bar when the 3*d 
bar is broken 

5 Figure 25 depicts a torque-time plot of healthy machine and one rotor bar-broken machine. 

Figure 26 depicts rotor velocities of healthy and shorted machines. 
Figure 27 depicts rotor torques of healthy and shorted machines. 
Figure 28 depicts rotor bar currents of shorted machine. 

Figure 29 depicts Kim and Bryant's bond graph of an induction motor with state variables. 
10 Figure 30 depicts angular position and velocity of rotor axis. 

Figure 31 depicts flux in rotor a windings ; the P winding flux is similar. 

Figure 32 depicts flux in stator a windings ; the winding flux is similar. 

Figure 33 depicts rotor velocity of a machine with a broken rotor bar. 

Figure 34 depicts stator current in the Frequency domain with broken bars. 
2 5 Figure 35 depicts torque-speed characteristics of the ideal and degraded machines. 

Figure 36 depicts power spectrum of the machine response and noise- 
Figure 37 depicts noise in the signal of the angular velocity of me degraded rnachine. 

Figure 38 depicts channel capacities with a broken bar. 

Figure 39 depicts rotor velocity of ideal and shorted rmchines. 
20 Figure 40 depicts power spectrum of angular velocity for the shorted machine. 

Figure 41 depicts spectral content of stator current of phase A; (a) Ideal naachine. (b) Shorted 
machine, (c) Ideal rnachine of [15]. (d) Shorted machine of [15] 

Figure 42 depicts spectral content of stator current of phase A with two severely shorted coils. 

(At - Rs2 =0.0079 a "Wi- nsi =10) 
25 Figure 43 depicts channel capacities with one shorted cofl. 

Figure 44 depicts channel capacities with two shorted coils. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Claude Shannon formulated a mathematical theory of communication. His groundbxeaking 
approach introduced a simple abstraction, me cornmunicarion channel consisting of a sender (a source 
of information), a transmission medium (with noise and distortion), and a receiver (whose goal is to 
5 reconstruct the sender's messages), see figure 1 . 

The transmitter injects messages from an information source into the channel. The receiver accepts 
a signal from the channel that contains the transmitted signal altered by the dynamics of the channel, 
and corrupted by noise added by the channel. 

An analogy is mads between a machine component or a system and a communications channel. 

10 During operation, information is sent as a signal over a communications channel from transmitter to 
receiver (See Figure 2). The signal over the channel is altered by limited dynamic bandwidth, 
nonlinearities and noise. The goal is for the receiver to extract and reproduce the message, despite 
distortions and noise. Design of communications systems is aided by powerful theorems of Shannon 
( 1949), which establish mmimum signal to noise ratios for error free transmission. 

1 5 A machine component (or system) accepts a "signal" from an upstream component, by its function 

alters mat signal, and men passes the "signal" on to the next downstream component In the analogy of 
this article, a machine is a communications channel. When operating properly, the "signal" from an 
upstream component is "received" by a downstream component Faults in the machine that disrupt 
functionality alter the "signal". Faults will be viewed as agents that alter system parameters or 

20 contaminate the signal with "noise". Unless the signal to noise ratio is kept sufficiently high, 
downstream components cannot "resolve" the "signal message" error free, and the machine 
malfunctions. 

. In performing a function, a machine,, component, or. system accepts a stimulus "signal" from 
another upstream component, alters that signal via its mechanical* fosction, and then passes the signal 
25 on to the next component. The signal contains information, which can be envisioned as a "message" to 
other components in the machine. The "message" relates to the function or intended operation of the 
machine or machine cornpouents. The mechanical fuiuiuon often includes kinematics of morion and 
dynamics of operation. 

Here we will strike an analogy between a machine component and a communications channel. The 
30 transmitter, an upstream component, activates our machine component "channel" with input signal x(t). 
v>ssage of the "signal" through the channel is associated with component functionality: compoiient 
kinematics and dynamics alter the signal. The component response defines the outpul y(t). When the 
component channel operates properly, the "message" contained in the signal "received" by downstream 
components can be unambiguously "resolved". Defects and degradation of the component afflict normal 
35 operation, "distorting" the signal and contanuhating it with "noise" n(t). Unlike electromagnetic 
communication channels, the signal may pass through multiple power domains: electrical, mechanical, 
solid, fluid, chemical, biological, etc. along its path through a component or a machine system. We can 
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view a machine as a channel consisting of individual component channels connected together to form a 
larger channel. 

The theory depicted and described in this application may be adapted in unique ways not 
contemplated by others, including, but not limited to, Shannon's theory. A component is designed to 
have functionality, which can be defined in terms of the (designer's) intended reaction of the component 
to an excitation. Degradation alters the component response. For example, fretting corrosion of the 
surface of an electrical contact changes the electrical impedance through the contact Although this 
alters the response to a voltage stimulus, the resulting signal distortions caused by changes of the 
electric contact impedance are often posed hi terms of an effective noise riding on the transmitted 
signal. Thus degradation of the contact via fretting is often modeled as an effective noise source and/or 
an impedance change. 

In conmiunications theory, Shannon's theorems traditionally estimate the ma x imu m rate of 
information C that can be transmitted through a communications channel, given its bandwidth w and 
ratio of signal to noise powers S/N. Designers of traditional comrnnnications channels considered C to 

15 be fixed, and their designerToxts focused on designing transmission or encoding schemes that would 
increase the rate of information R up to its upper limit, the channel capacity C. If applied in a 
nontraditional manner to ma chiner y, Shannon's theorems can yield a threshold signal to noise ratio 
(S/N)t. In the commumcations channel analogy, dynamics inherent in the component functionality can 
be included in bandwidth w and in the signal to noise ratio S/N of the channel capacity C. These 

20 dynamics may change as the component degrades, causing C to change. For a typical machine and 
components, signal transmission rateR should be constant, since machine or component operation is 
often repetitive (or periodic) and at or near steady state: the machine controller and/oi upstream 
components continue to inject their si gnals into a machine (or component), regardless of its condition. 
When applying communications thsory to a machine component ox system, we will first trace the 

25 path of signal power flowing through a healthy (functional) component or machine, to define the 
communications channel through the component or machine. Along the signal path we will list the 
various forms of energy oi yo*er into which the signal is transformed Functionality will be defined in 
terms of the input to output response, for components or the system: if for a given set of input 
excitations, the output response matches within some tolerance the desired output, the component or the 

30 system is functional; otherwise it is dysfunctional. If needed, we will consider each separate 
energy/power domain and its transduction as a commumcations channel, and then connect these 
channels together in a manner consistent with the machine's runctionality and design. Bond graphs 
(Kamopp, Margolis, and Rosenberg, 1990), which map power flows through dynamic systems, can be 
useful, since bond graphs readily handle systems with diverse energy domains in an energetically 

35 consistent maimer. 

After analyzing the healthy system, we will then incorporate component faults and degradation 
modes into the system model. To affect functionality, the degradation effects must alter or block the 
flow of signal power through the component. Questions we must answer include: How does each 
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degradation mode alter the signal flow, and affect system or component parameters? Does the particular 

degradation cause components to become nonlinear? Does the particular degradation generate another 

signal, Le., noise? We will incorporate degradation into the system model as changes to existing system 

model parameters or as additional elements (e.g., sources of noise). Location of each degradation mode 

5 in the system model will be consistent with the locus of the degradation in the physical machine 

component or system. 

Aspects of this method include: 1) Individual components, or an entire machine system consisting 
of multiple components, can be analyzed. 2) System malfunctions can be predicted, including individual 
faults and those due to a collection of seemingly healthy components. Enois from slightly degraded but 
10 individually healthy components can stack through a machine system, rendering it unable to meet 
tolerance. 3) The current status of the system, and time to system malfunction can be estimated by 
simulations based on these models. 

Other aspects of the diagnostic procedure include: Determine and trace the path of tie signal flow 
through the healthy system, from signal in to.signal out. For the sick system, model the faults with noise 
15 sources or parameter changes. Multiple system outputs may exist At each output, tally the signal power 
and the total noise power to obtain a signal to noise ratio S/N. Estimate the bandwidth w for the signal 
path through the degraded component coinmunications channel, using the enhanced system model. 
Apply Shannon's theorems to diagnose the absolute health of the machine component or machine 
system The health of each individual component in a rrachine system can be assessed, and likewise the 
20 health of the entire machine system. 

The analysis of each machine component commimications channel may contain the following: 
. worthy Machina m*Hp! which has no faults and functions perfectly. This is an idealization that 
reflects the machine designer's original intended concept The output y 0 (t) of a signal x(t) 
propagating through this ideal machine will define the intended message or signal y 0 (t) that the 
25 machine component or machine system corrmnmications channel is supposed to transrnrt and 

receive. The signal powers S - P{y}and Si - P{y 0 } defined by S = P{y}= 'flfeqn (1), 

both in the channel capacity C and the rate of iriformation transmission R through the channel, will 
be based on this ideal machine. The resulting model is simple, the concept of perfect health is well 
defined, and the signal y 0 (t) that the receiver is supposed to receive is well defined. 

30 

- Machine Emits These include common degradation faults for a given component Common 
examples include pitting of gear teeth, fatigue cracking of shafts, and deterioration of insulation on 
electric motor stator or rotor coils. 

35 . Madura Fault Models incorporate the Machine Faults as sources of noise n(t) and/or changes 
in system parameters consistent with irnr^rfections, faults and degradation modes of a particular 
machine element Noise will be defined as any signal component that should not be in the perfectly 
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transmitted and received message signal y 0 (t). This may includes haimonics generated by nonlinear 

elements. 

. ng qrariPri M ar tyr" Model : This is the overall system model that results from adding the 
5 Machine Fault Models to the Healthy Machine Model It includes sources of noise n(t) and changes 

in system parameters. When all noise sources are zero, the healthy machine results. Transmission 
of the signal x(t) through the degraded machine (noisy cominunications channel) induces received 
signal y*(t), generated by signal x(t) (sent through as y 0 (t))and noise n(t). 

10 The analogy may also be extended to a set of machines, a process, a manufacturing or assembly 

method, or others. The analogy may hold for a series of "information channels" as seen in figure 3. 

The model-based diagnostics is based on fundamental first principles of physics and 
information theory. The methods uses sensor signals to tune the parameters of a model of the system, 
such that the model then mimics the operation of the real system. Diagnostics are performed on the 

15 model. The diagnostic system can be designed as part of the design of a new machine. Also, models 

allow what if predictive scenarios for a machine's possible furore health and functional condition,' given 
ce rtain observed trends in the machine's health. For a different machine or a new design, only the 
model of the operation of the system must be altered, not the basic diagnostic algorithm. Models also 
avoids interpretation of complex sensor signals, trying to figure out what a particular peak or dip, or a 

20 band of frequencies in a signal means, in terms of machine health. Instead, time wise changes to 
machine parameters can be followed, and projection of these trends can be used to forecast future 
health. Models also permit incorporation of knowledge of faults, and the intent of the designers of 
machinery, into the diagnostics routines. 

To quantitatively analyze transmission through the channel, "Shannon introduced a measure of 

25 the amount of information in a message. The measure is related to the probability of occurrence of the 
events for which the messages are about A message that informs the receiver that a rarely occurring 
event is about to ^ppen contains the most information. A message informing about an already "known" 
event conveys little information. Information entropy, a measure of the average amount of information 
(or uncertainty) in a message, can be defined [1] as 

n 

30 h=\ j Bl ' 1 o) 

- x poftto&i p&yt* % xis c ° ntinuous 

Here p t is the probability of occurrence of the message's event xi if the random variable is discrete, and 
p(x) is the probability density function for the random variable * , if the random variable x is 
continuous. Here Pi is the probability of occurrence of the message's event *,xi if the random variable 
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is discrete, and />(*)is the probability density function for the random variable * , if the random 

variable x is continuous. 

Shannon's entropy rate (R) measured a source's information production rate, and the channel 
capacity (Q measured the information carrying capacity of the channel. As per one of Shannon's 
5 theorems [1], if R ^ C , then there exists a coding technique which enables transmission over the 
channel with an arbitrarily small frequency of errors. Ims restriction holds even with bounds the noise 
in the channel. A converse to this theorem states that it is not possible to transmit messages without 
errors if R > C . Thus the channel capacity is defined as the maximum rate of reliable information 
transmission through the channel. 

In another theorem, Shannon derived the channel capacity for a time continuous channel with 
additive white Gaussian noise. His expression 



10 



(2) 



involves the average transmitter power, 

1 0 



15 of the signals x o (*) , the power of the noise, 
1 o 

and bandwidth & of the channel in hertz. If the bandwidth is non-flat, then the capacity of the channel 
is given by 



df. (5) 



20 S imilar ly, the entropy or information rate for messages 

r ^ < 6) 

derived by Shannon involves S t , the average power of the desired signal to be transmitted, N g , the 
maximum allowed RMS error between recovered and original messages, and CO i , the signal bandwidth. 
Shannon's communication theory could be applied to the fault diagnosis of machine systems. 
25 A machine component (or system) accepts a signal from an upstream component, by its function alters 
that signal, and then passes the signal on to the next downstream component In Bryant's analogy, a 
machine conveys information in a signal and is thus a communicanons channel. When operating 
properly, the signal passes through the system and is successfully received within desired tolerances at 
the machine's output. Faults that disrupt operation alter the flow of signaL Faults will be viewed as 
agents that contarninate the machine's signal with "noise". Unless the signal to noise ratio 



30 
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kept sufficiently high, downstream components cannot resolve the signal message error free, and the 

machine malfunctions. 

Noise is denned as an "unwanted signal tending to obscure or interfere with a desired signal", 
as "any signal which interferes with the transmission of a signal through a network or tends to mask the 

5 desired signal at the output terminals of the network' 1 , and as "an unwanted signal tending to interfere 
with a required signal". Thus noise is the difference between the actual signal received, and the signal 
desired to be received To apply this definition to mechanical systems, we must define the desired 
signal. We shall call mis desired signal the "ideal" signal X 0 (t) . Note that x o (f ) , an idealization, must 
be produced by a system without noise. This is possible only with models, not with real systems. 

1 0 The "ideal" and "degraded" models may be defined as follows: 

. The ideal machine model has no faults and functions perfectly. Its output defines the signal x o (0 that 
the machine channel is supposed to receive. From this, we can estimate signal power S { , 
• The degraded machine model is the overall system model mat results from adding faults to the model. 
We will incorporate feints as noise »(0 . Thus the signal *(0 = * 0 (0 + "(0 contains noise n & ■ 

15 denned as any signal component that should not be in the perfectly received message signal. Noise is 
any deviation from the ideal signal, including unwanted harmonics generated by nonlinear elements. 
This will estimate the noise power N . 

We can incorporate degradation or imperfections into the system model Degradation can be instilled 
in a bond graph model by varying bond graph parameters, adding noise (effort or flow) sources, or 

20 changing the power pathways. 

The models may take various forms. These forms may be any form appropriate for use in the system 
"of application. For example, these forms may be heuristic, neural netwcm^rdetennmistic, probabilistic, 
and others. 

The method and model system may be implemented on a computer system, S see Figure 4). 

25 The term "computer system" as used herein generally describes the hardware and software components 
that in combination allow the execution of computer programs. The computer programs may be 
implemented in software, hardware, or a combination of software and hardware. A computer system's 
hardware generally includes a processor, memory media, and input/output (I/O) devices. As used 
herein, the term "processor" generally describes the logic circuitry that responds to and processes the 

30 basic instructions that operate a computer system. The term "memory medium" includes an installation 
medium, e.g., a CD-ROM, floppy disks; a volatile computer system memory such as DRAM, SRAM, 
EDO RAM, Rambus RAM, etc.; or a non-volatile memory such as optical storage or a magnetic 
medium, e.g., a hard drive. The term "memory" is used synonymously with "memory medium" herein. 
The memory medium may comprise other types of memory or combinations thereof. In addition, the 

3 5 memory medium may be located in a first computer in which the programs are executed, or may be 
located in a second computer that connects to the first computer over a network. In the latter instance, 
the second computer provides the program instructions to the first computer for execution. In addition, 



013270.00015:124601.01 



WO 02/27418 



PCT/US01/29935 



~ 14 - 
the computer system may take various forms, including a personal computer system, mainframe 

computer system, workstation, network appliance, Internet appliance, personal digital assistant (PDA), 

television system or other device. In general, the term "computer system" can be broadly defined to 

encompass any device having a processor that executes instructions from a memory medium. 

5 The memory medium preferably stores a software program or programs for the reception- 

storage, analysis, and transmittal of information produced by an Analyte Detection Device (ADD). The 
software program(s) may be implemented in any of various ways, including procedure-based 
techniques, component-based techniques, and/or object-oriented techniques, among others. For 
example, the software program may be implemented using ActiveX controls, C-H- objects, 7avaBeans, 

10 Microsoft Foundation Classes (MFC), or other technologies or methodologies, as desired. A CPU, such 
as the host CPU, for executing code and data from the memory medium includes a means for creating 
and executing the software program or programs according to the methods, flowcharts, and/or block 
diagrams described below. 

A computer system's software generally includes at least one operating system such Windows 

1 5 NT available from Microsoft Corporation, a specialized software program that manages and provides 
services to other software programs on the computer system. Software may also include one or more 
programs to perform various tasks on the computer system and various forms of data to be used by the 
operating system or other programs on the computer system. The data may include but is not limited to 
databases, text files, and graphics files. A computer system's software generally is stored in non-volatile 

20 memory or on an installation medium A program may be copied into a volatile memory when rumiing 
on the computer system. Data may be read into volatile memory as the data is required by a program. 

Further, the method may be implemented across a set of networked devices (See Figure 5). 
The method may be performed remotely from the system Further, the results of the method may be 
transmitted, stored, processed, and accessed across a network, among others. 

25 For example, parameters for a model of a patients health may be stored on a smart card. 

These may be accessed arH^nmbined with the method to determine a change in state of the patient's 
health. In another exemplary embodiment, a machine may be located in a remote location. A service 
provider may periodically access data from the machine from a remote location and diagnose the 
machine. These diagnoses may be used in predicting the failure of the machine. Further, these 

30 diagnoses may be used in placing an order for a replacement. 

Figure 6 depicts a flowchart for diagnosing according to the invention. The method may be 
implemented in software and/or hardware. Further the method may include some or all of the steps in 
various combinations. 

In a first step, the user is directed to assemble detailed dynamic systems models of the rnacbine 
35 system in question. The models may possess a one to one correspondence between physical components 
or elements in the real physical system, and elements in the dynamic systems model. One may include 
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all possible faults and potential failures in the system models. This invention may tabulate the effects of 
feulte as "noise" in the system. Noise in a signal is the difference between the actual signal and the 
expected signal. In the model, noise may be induced by changes in parameters of dynamic system 
elements, which men alters any signals passing through a system Or, if a certain fault cannot be 
5 described by these means, then sources of noise (often white noise) will be inserted into the system 
model, at locations in the model that are consistent with the locations of the faults in the real machine. 
The intensity of these noise sources can then be adjusted to make the model behave like the real 
machine. 

One may then judiciously monitor the machine or system Excite the machine or system, and 
10 observe the machine's or system's response, for example, via the sensor outputs. 

One may then tune the model's parameters, so the model mimics me real system Excite the 
system model with the same excitation as the previous list item Compare the outputs of the model to 
the corresponding outputs of the real machine or system If the model's outputs differ from the real 
machine's or system's outputs, adjust or change values of model parameters, including the intensity of 
15 the noise sources, until the model's outputs closely match the actual system's outputs. 

One may men manipulate the model, which now mimics the real machine or system in its 
present condition: 

From the model, one may calculate the channel capacity, C, of the machine. C is the inaximum 
amount of mforrnation that can be observed successfully conveyed through the machine. The channel 
20 capacity depends on the design and construction of me system, and the present condition of the system, 
which results from manufacture, aging and damage. Faults are encoded as "noise" in the model. 
Analytically, C depends on the strength of the noise levels in the system, relative to the strength of the 
excitation system response signal. 

For a desired job to be performed by the machine, one may calculate the rate of mforrnation R 
25 associated with the job .R depends of the desired speed at which the maclime does the job, the desired 
loads, the complexity of the job, and the desired accuracy at which the rriachine should do the job. R is 
measured in bits of information per second. 

- Compare R to C. If R < C, the machine will perform the desired job within the desired 
precision If not, the system has runctionally railed. 
30 The comparison of R to C may yield a diagnosis. Alternately parameters of the tuned model 

may yield a diagnosis. Further, this diagnosis may be associated with the determined noise. In 
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addition, the noise and/or diagnosis may be indicative of combined faults. Further, combined variances 

in parts, while within tolerance limits, may comprise a fault, defect, or otters. 

The method may be repeated over time to build a prognosis of the machine or system For 
example, a prognosis may predict the failure of a part 
5 Further, the method may be applied to many systems such as those depicted above. In 

addition, parameters from the tuned model may indicate the type or state of a defect, fault illness, or 

condition, among others. 

In typical applications, the method may involve formulation of extremely detailed models of 
machine devices to describe a machine's condition. These are critical to success. For example, included 

10 are bond graph based models of a motor, a gear box, and other mechanical transmission components. 
These extremely detailed models: (a) exhibit a one to one correspondence between elements in the 
model and components in the real system; (b) incorporate all known effects of the device into the 
model, including defects; (c) include in the models via finite element concepts instilled into bond 
graphs the dynamically distributed nature of components in the real system, and (d) use noise sources to 

15 account for defects and degradation of components. Simulation of the motor and gear box models can 
generate the complex spectra measured during operation of these devices. 

The models tabulate the effects of system faults (machine maladies) in a very novel way: as 
'noise" in the machine. Noise is the difference between the actual signal received, and the expected 
signal that should be received. As a rmchine degrades or ages, the difference between actual and 

20 expected signals becomes larger, and thus floise levels increase. These noise methods permit 

incorporation of faults into the models that heretofore could not be described analytically. The concept 
of noise has been used heavily in electronics and communications engineering, to design around noise 
"faults" always present in these electronic and electromagnetic systems. Electronic noise, including 
resistor noise, shot noise, burst noise, and flicker noise among others has been generally tabulated or 

25 modeled with noise sources placed in a model of the electronic circuit This work imported this body of 
knowledge to mechanical, hydraulic, and other physical systems, but in addition, systems extended the 
modeling scliemes of noise to include noise induced by changes in parameters of the system. 

The method also applies techniques of information theory to machinery - as opposed to present 
applications that are limited to electronic communications systems - to quantitatively assess the current 

30 health state of a machine. The method treats a macinne, such as a CNC engine lathe, as a noisy 
corrmurnications channel, to assess reliability and functional condition. A message transmitted and 
received over a cornmunications channel picks up noise due to imperfections present in the physical 
channel. For example, music transmitted over an AM channel is overwhelmed by buzzing when the 
receiver is near electrical power transmission lines: the transmitted musical message is obscured at the 

35 receiver by electrical noise. In an analogous manner, a machine transmits a message over a machine 
channel. For example, a lathe, viewed as a communications system, has transrnitter = CNC controller, 
channel = (drive motor + gear box + lead screw, + tool carriage on ways + aitting tool / workpiece), 
and receiver = workpiece. "Noise" includes effects of fatigue, spurious vibration (from other machines), 
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and other errors due to wear of machine and cutting components. A transmitted "message" is properly 

"received" if the finishe d part is within tolerance, or in a general machine, if the machine performs its 

function within specified tolerances. Excessive noise in the machine systemmay cause a part to be out 

of tolerance, or causes the general machine to operate outside the specified tolerance limits. With this 

5 view, Shannon's comrrmriications theorems may be applied to machinery'. Shannon's theorems may 

accurately estimate the limits on the amount of information per unit time C that can be sent through a 

noisy communications channel. C depends on the channel's state, including dynamics and signal to 

noise strengths (ratios). For a lathe, making a part of certain geometric complexity at a given speed, to 

within a desired (fidelity) tolerance is characterized by an information rate R. JfR < C, Shannon's 

10 theorems predict success; if R > C 5 the part will be out of tolerance. As a machine deteriorates, C 
decreases, and eventually R > C. Now the machine cannot make the part with the same speed and 
tolerance. The channel concept appears to be a very sensitive discriminator of a machine system, even 
for the stacked effects of a collection of moderately degraded components. 

The method may be used for predicting the future conditions of machinery, for scheduling 

15 maintenance and avoiding functional and catastrophic failures of said machinery. The method can 

forecast if a complex system is capable of doing a given task, within a specified tolerance. A multitude 
of parameters associated with the machine's model may be tuned, such that the model emulates the real 



These modeling and system assessment techniques could be useful to designers of machinery, 

20 to assess the efficacy, reliability and durability of a design under various user conditions. 

In addition to mechanical systems, these methods could apply to almost any kind of dynamic 
system, including chemical, electrical, medical, manufacturing and processing, and business operations, 
among others. For example, in the medical world, a detailed model describing the dynamics of the 
caidio-vascular system could be developed. This model would possess multiple parameters that 

25 describe behavior and condition of the heart and blood vessels, and their interactions with other body 
systems such as lungs and kidneys. The model could be tuned from medical signals and data derived 
from tests and procedures, such as Electro-Cardiogram, blood pressure, and data from ku tests and 
radiology. After timing the models, a channel capacity C could be estimated to assess the condition of 
that system, and compared to a rate of inforrnation R. This comparison would assess the health state of 

30 the patient The rate of information would describe the ability of the cardio-vascular system to perform 
at various levels characterized by task speed, load, complexity, and tolerance. Since the rate contains 
these factors, degrees of health and sickness could be assessed quantitatively or assessed, in a formal 
manner. This could automate medical diagnostics. Medical prognostics would extrapolate trends of 
parameters in the model, or trends contained in the data, and apply the channel capacity and rate of 

35 information concepts of communications theory, to forecast future health scenarios. 

These methodologies could be extended to evaluate business practices, procedures, and 
enterprise structures. A business operation has dynamics imposed by its processes, people, and 
structure. The application would treat an enterprise as an imperfect communications channel, and 
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construct models of information flow through mat system Transmitters— the orders— will send 
information over imperfect "enterprise communication channels". Imperfecnons—problems in the 
enterprise, or interference between conflicting missions-adds "noise" to channels. Receivers— the 
customers— must receive the message— the product— within tolerances— customer expectations— 
5 despite noise. The application would define "channels" through enterprise units, construct models that 
mimic these channels, and then apply communications theory to diagnose and prognose these channels. 

The models in these embodiments and claims can take various forms: from structured methods 
such as bond graphs, differential equations, and finite elements, among others, to heuristic methods 
such as neural networks, fuzzy logic, expert systems, and other computer methods. 

10 The method for applying communication theory to machines and systems need not be limited 

to signals derived from models. The method could be extended to signals measured from real systems. 
Here the ideal signal Xo(t) could be approximated from measurements taken from a real machine, or 
from several rnachines, in excellent condition. The difference between Xo(t) and the signal x(t) measured 
from a degraded machine could replace those derived from models, mentioned earlier. Similarly, the 

15 difference could be used to confirm that a machine operates within tolerances. Further, an ideal signal 
could be a signal from a machine with a known defect. The difference between the signals would then 
confirm a specific defect, among others, 

pypmplary applica tion tn a squirrel caqft indlictinn motor 

20 Equation numbers in this example refer to equations listed in this subsection. Similarly, an 

appendix is attached that is referenced in mis subsection. 

One exemplary application of the invention is a method for diagnosing an induction motor. 
For example, a motor has two major sub systems: a rotating rotor and a static stator. Induction machines 
can have a wound rotor, or a squirrel cage rotor. Widely used squirrel cage induction machines exhibit 

25 great utility for variable speed systems and are simple, rugged, and inexpensive. The squirrel cage rotor 
is a structure of steel core laminations mounted on a shaft, with solid bars of conducting material in the 
rbtoYslots, end rings, and usually a fan. In large machines, the rotor bars may be of copper alloy, dri vcu 
into the slots and brazed to the end rings. Rotors of up to 50 cm diameter usually have die-cast 
aluminum bars. The core laminations for such rotors are stacked in a mold, which is then filled with 

30 molten alumimim. In this single economical process, the rotor bars, end rings and cooling fan blades are 
cast at the same time. 

Figure 7 is a schematic of a squirrel cage induction motor. A substantial literature modeling 
induction motors employs Park's (1929) two-reaction theory, which accounts for magneto-mechanical 
energy transduction via multi-port inductances. From Park's model, Ghosh and Bhadra (1993) 
35 formulated the bond graph in Figure 8. We altered Ghosh and Bhadra's bond graph to partition and 
make explicit the electrical, magnetic, and mechanical energy domains; to form a one to one 
correspondence between physical components in the machine, and elements in the bond graph; and to 
append additional elements to the bond graph to make it more consistent with real induction motors. 
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When energized by an AC supply voltage, the stator coils form a radial magnetic field vector 
that rotates within the interior of the stator, about its central axis. Within this interior the stator field 
cuts through the squirrel cage rotor, including conductor bars mat extend airily. This time varying 
field induces a voltage over the rotor bars. Resulting bar currents flow in the sequence: bar -» end ring 
_> opposite side bar opposite end ring original bar. Induced by this time varying current loop is a 
secondary magnetic field, which attempts to align with the stator field. However, because the rotating 
stator field induced the secondary field of the rotor, the stator field leads the rotor field, and 
consequently, the rotor chases the stator field, always following. This is motor action (Xawrie, 1987). 
Tne induction motor speed depends on the speed of the rotating stator field. 

The real system we will consider is a two pole, 'Y' connected three phase squirrel cage 
induction motor. In (Ghosh and Bhadra, 1993; Sahm, 1979; and Hancock, 1974), a multi phase 
induction motor was modeled with an equivalent two-axis representation. Each phase winding generates 
its own magnetic field, which can be represented as a vector aligned along me axis of the wmding. The 
sum of these phase vectors produces a phasor vector. If the phase vectors vary properly with time, the 
phasor rotates. 

A transformation from three phases (a,b,c) to two phases (a,p) was represented in (Hancock, 
1974) in matrix form. If the 'a ; and 'a' phase wmdings are co-axial, the induced Magneto Motive 
Forces (MMF) of the V and V phases of the three and two phase systems are co-directional By 
appropriate changes to the two phase currents, the magnitude of the phasors of the three and two phase 
systems can be made equal. Ghosh and Bhadra (1993) represented this in their bond graph via 
reformer elements in the stator section. The two phase currents were represented in terms of three 
phases as 

[V2/V3 -lA/6 -l/V6~ 

o i/>/2 -i/ir 
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Under assunqjtions of a spatially sinusoidal disliibution of MMFs, and ignoring magnetic losses and 
saturation, Ghosh and Bhadra (1993) expressed a symmetric mdncTRm motor in an orthogonal 
stationary reference fiame with a and p phases fixed on the stator as 
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Equation (2) relates stator voltages'to stator and rotor currents. In addition, needed is the electro- 
magnetic motor torque for a P-pole machine, expressed as 



This motor torque is balanced against other torques via 



(3) 
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T.'J&f+co.+Z (4) 

Terms on the right side of equation (4) represent rotor inertial torque, shafVbearing damping 
torque, and load torque, respectively. In equations (2) to (4), and V p , are a and p axis stater 
voltages; and ip s are a aad 0 axis stator currents; and i pr arc a and (} axis rotor currents; R s and R r 
5 are stator and rotor resistances; L,, Lm and Lr are stator self inductance, mutual inductance and rotor self 
inductance; T e and Tl are electro-magnetic torque and mechanical load torque; J is the moment of 
inertia of the rotor, c is the viscous resistance coefficient; Or and com are electrical and mechanical 
angular velocities of the rotor; and P is number of pole pairs. 

Ghosh and Bhadra (1993) represented equations (1) to (4) in their bond graph, reproduced in 
1 0 Figure 8. They used modulated gyrators MGYir^Lmip,, MGYr^LriBr, MGYirs^Lmioa* MGYx4=Lrior, 
to represent the electro-magnetic torque of equation (3); employed transformers TF:mi, TFnrfe, TFrmj, 

TF:m4, TF:m5 with moduli • n\ = ^j, m 2 = m z = — n/6, m 4 = ~j2,m 5 = — to implement the 

mathematical transform of equation (1); and excited the system with effort sources MSe:v a> MServb, and 
MS c :v c having sinusoidal voltages with equal amplitudes but 0, j&G, and 2n/3 phase lags, respectively. 
1 5 Although mis correctly programs the governing equations for a three phase induction motor, it lacks a 
correspondence between bond graph elements and real system components. Moreover, elements and 
theiT constitutive laws involve only electrical and mechanical energy domains. Faults or design 
parameters relevant to the magnetic domain are only implicit in the mutual inductances, posed as 2 port 
inertances La and I:p with constitutive laws 
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(5) 



In Xc», 3^3, Xefy and are flux linkage of the respective windings. In Figure 8^ five integral 
(independent) causalities exist on inertance energy storage elements, with system state variables 
Aps, Xor, V» m ^ k wnere h is the rotor angular momentum 

To represent real system elements or components explicidy, certain bond graph elements 

25 should be moved, altered or added. In Figure 8, a and (3 phase stator resistance elements, Rsa and R S|J 
should be split into three stator coil resistances R^, R*b, and R«, without alterating the governing 
equations. The revised bond graph shown in Figure 9 moved R^ and R«p back through the transformers 
in front of the phases. To maintain an equivalence between Figure 8 and Figure 9, we must relate Rsa, 
Rib and R a to R* 0 and R«p. Since most motors possess symmetry between phases, let Ka = R$b = R sc = R, 

30 and = Rs0 = R s . For the bond graphs of Figure 8 and Figure 9 to be equivalent, the voltages (efforts) 
to the 2-port inertances on the stator sides must be equal in both Figure 8 and Figure 9. The causality in 
both Figure 8 and Figure 9 asserts that these voltages to the 2-port inertances arise from the neighboring 
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Inunctions. Summing voltages from other bonds to these l-junctions, and equating these respective 
voltages between Figure 8 and Figure 9 gives 
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5 By solving for i^, we obtain equations in temis of resistances and transformer moduli 
i fi (mjn t ml + m^m^R 



(8) 



By replacing the transformer moduli, mr-n* of the three phase to two phase transformation with real 
numbers, n h = ^.m, = m> = m 4 = V5, m 3 = -^2 , which is given in equation (1), we find that 
R,=R, i.e., Rto= Rsp 3 1^ = ^ = R 3 c 

Simplified Representation of the Signal and Modulated GY Elements 



In terms of the 2-port I field of equation (5), equation (3) can be rewritten as 



(9) 



From this relation, Figure 8 can be rearranged into the form of Figure 10, where the modulated gyrators 
MGY: r« = and MGY: x l2 - are modulated by the flux linkages \* and V of the 2-port 
inertances. 

The number of squirrel cage rotor bars depends on the rotor's size, and usually, tens of bars 
are in one rotor. In this study we consider the squirrel cage rotor with five bars (numbered 1 to 5) 
depicted in Figure 1 1. Shown also is tie rotor magnetic field (dashed line), with north poles (N) on top 
of the rotor, and south poles (S) beneath, and bar currents. Currents directed out of plane are denoted by 
a V, and currents flowing into the plane arc denoted by a V. Each end of each rotor bar is attached to 
a solid end ring. Induced currents flow through each bar and end rings. With five bars, there exist five 
different currents (flows) in this rotor. At the instant of the rotor position shown in figure 1 1 -(a), the 
sums of the currents induced by the rotating magnetic field of the stator in bar 1, 2 and 5 must be equal 
to the sum of the currents in bar 3 and 4. Likewise, the current sunimation of bar 1 and 5 at the position 
of Figure 1 1 -(b) must equal the sum of currents in bars 2, 3 and 4. In Figure 1 1, the thickness of each x 
and • shows the relative current magnitude in each bar. 
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To incorporate individual rotor bars into the bond graph, the a and p phase currents and 

voltages of the rotor should be split into separate bar currents and voltages. The a, b, c and cc p axes 

are stationary with respect to the stator, but because the rotor rotates relative to these axes, bar currents 

must depend on the rotation angle 6 of the rotor. Using results in Hancock (1974), rotor bar currents 

can be related to the a, (5 phase currents as 



(10) 



In equation (10), i* represents the current in the k* rotor bar (k = 1, 2, ... n), 7^ and X 3 r are rotor 
currents from Figure 8, and magnitude modulus m depends on the total number of bars, n. For n = 5 
bars, we will have currents in to i*. Accordingly, rotor bars can be incorporated into the bond graph of 
Figure 10 via 9-modulated transformers. 

Figure 12 shows the tra^formation of a and p phase currents into individual rotor bar 
currents, where the transformer rnoduli are 



mr k = mcos<0 + > 



k = 1.2,....n 



with n = 5 . 



(ID 



(12) 



In Figure 12, the battery of 0-junctions on the right side completes the summation of a and p phase 
currents demanded by the right side of equation (10). The voltages that sum over the two 1 junctions 
located between the I fields and the MTFs give rise to 

K + (nwjH + ( mr 2 )^ + * * ' + ( mr s)^ = 0 
V + (mr 6 )Z 6 + {mr^ + • • • + (mr I0 H o =0 

Here the flux linkage \ u X2, . . Xu associated with rotor bars are located to the right of the MTF's. To 
obtain the torque contributed by each bar, equation (10) for k - 1, 2, .. ., 5 is rewritten in matrix form as 



(13) 



VI 
V5_ 



cos& 



sinO 



cos( 0 + 2l) S m(0+y) 
cos(0 + — ) sin(0 + — ) 

cos(0+— ) sm((9^— ) 



Le. 



irotor^Aifwo phase 



(14) 



The two column vectors of the 5x2 transformation matrix A form an orthogonal set for any value of 
rotor rotation angle Q; the rank of A is 2. For the mxn (mOn) matrix A having rank n, there exists 
(Strang, 1 988) an nxm left-inverse B such that BA=I n , where I n is the identity matrix of order n. In our 
model 



A T A = 7W 2 — 
2 



0 



(15) 
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and the left-inverse of A is A T if m 2 =|, i.e., the transformer modulus m has a value which 

normalizes A T A. For a rotor of n bars, m = ^. The proof is shown mAppendix for tms subsection. 
From equations (14) and (15), the inverse transformation is 



V2 

'Is 



COS* cos{eA coswA costfA costfA 

8 

srnff sin(0-A sin(0+^) sin(0+^) sra(0 + -f) 
5 5 3 J J 



(16) 



10 



15 



20 



25 



If substituted into the rotor output torque equation (9), the electromagnetic torque becomes 

2(*-l)* 



5 p sj2 



^cosj 0 + 



(17) 



The revised bond graph in Figure 13 includes stator and rotor bar interactions based on equation (17). 
Here the moduli of the k* modulated gyrator is 

2(fr-l>r V 



(IS) 



where n - 5 for Figure 13. Finally, the electric resistances of the rotor were grouped with eachiotor bar 
in s manner similar to mat of the stator resistances. 

The bond graph in Figure 13 models the interaction between stator coils and rotor bars with 2- 
port I elements— inductances— in the electrical energy domain. An inductance only describes storage of 
magnetic energy. Neglected are power losses and leakage effects 'in the magnetic domain, which may be 
caused by component deterioration. To describe these interactions, we replace all I inductance elements 
with equivalent combinations of gyrators and C elements, without violating causality. Figure shows 
equivalent bond graph representations between an'Tand a GY and C combination; and a TF and GY 
combination. 

In Figure 14, n is the gyrator modulus (the effective number of coil turns); m is the transformer 
modulus; I is the flux linkage; (jj is the magnetic flux [Wb]; M is the magneto motive force [A]; p is 
the permeance of the magnetic circuit element [H]; ei and e 2 arc efforts; and fi and f 2 are flows. In 
Figure 14 -(a), through the gyrator relations A = nq> and W = M . Using the constimtive law of the C 
element, M =» $lp, the two port I elements pertaining to the a and p phases were converted into 2-port 
C elements that now represent interactions between magnetic flux and magnetomotive force of the 
stator and rotor. Figure 15 shows the new bond graph with five rotor bars and the GY - C - G Y 
combination that replaced the 2-port I. The gyrators were then moved through the bond graph to new 
locations more consistent wiih motor components. The GY to the left of the 2-port C was moved into 
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the electrical section, where it now represents the action and number of turns of the stator coils. The GY 
leap-frogged the transformers that were based on equation (1), changing moduli of these transformers 
according to Figure 14 -(b). The GY to the right of the 2-port C skipped over a 1-junction, converting 
that 1-junction into the 0-junction shown in Figure 15. Similarly, a 0- and 1-junction to the left of the 2- 
port C in Figure 13 were converted to a 1- and 0-junction in Figure 15. In the bond graph of Figure 15, 
electrical energy inputs, transformation of energy from electrical domain to magnetic domain, 
mathematical phase reformations, power interactions between stator and rotor bars in terms of 
magnetic flux and magneto motive force, and mechanical rotor output are all represented and labeled 
In Figure 15, the two sets of gyrator moduli n, and n. stand for the effective coil turns which relate 
electrical and magnetic variables of stator and rotor, respectively. 

State equations were derived from the bond graph of Figure 15 with n«i = n»2 = n& - n*. - 
R sb = Rsc - Rs- In terms of magnetic variables, the state equations are 



n s m, n s n^ 

TflfJ 



(19) 



m 6 J 



lb where the magnetic state variables are stator and rotor phase fluxes <|> a «, fo. <t>„r» and rotor angular 
momentum h. The constitutive law of the 2-port C element is 

(20) 



p, p m 
p p 



In the state equations, 

i 



i 



44 -l; n s \-p s P r -P m - 



(21) 



20 where the permeances P s - \,P m -^J are expressed in terms of coil turns and 

n, n s n r n r 

inductances of stator and rotor. Here n s is the n umb er of effective stator coil turns, n r the number of 
effective rotor coil turns, § the magnetic flux [\Vb], M the magnetomotive force [A], P is the Permeance 
[H], and h the angular momentum [N-m-s=kg-m 2 -sec]. 
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Simulations of asquinel cage mducnon motor used the bond graph shnulationtool, 20-SIM 
(Control Lab Products, 1998). For integration of state equations, a lUinge-Kutta 4 th order method was 
adopted. Values of the system parameters for the simulations are presented in 

Table 1 , some were identical to those used by Ghosh and Bhadra. 

5 Table 1 System parameters of a two pole, three phase squirrel cage induction motor 



10 



15 





Stator coil resistance 


0.078S 


Rri, IW..,R*io[Q] 


Rotor bar resistance 


0.0408 


V a , Vb, v c [V] 


Input voltage amplitude 


230 


[Hz] 


Input voltage frequency 


60 


L S [H] 


Stator inductance 


0.0153 


U[H] 


Rotor inductance 


0.0159 


L*[H] 


Mutual inductance 


0.0147 


n 


Number of rotor bars 


5 


ns 


Number of effective stator coil turns 


100 


Br 


Number of effective rotor coil turns (bar) 


1 


c[N-s/m] 


Mechanical resistance 


0.15 


Jfkg.m 2 ] 


Mechanical inertia 


0.4 



Shown in Figure 16 and Figure 17 are plots of rotor angular velocity and stator currents versus time. 
The rotor velocity rises slowly to a steady state value of about 377 rad/sec; the stator currents oscillate 
at the input frequency with initial large amphtude. After about 1.5 seconds, the motor reaches steady 
stare: the currents in stator windings decrease to a steady value and no osculation of rotor velocity 
exists. Figure 1 6 plots the rotor axis angular velocity vs tune when 230V, 60Hz three phase AC 
voltages are input to the stator coils. Theoretically, when 60Hz alternating inputs are given to a two 
pole AC motor, the output velocity should be 3600RPM (377 rad/sec) and the simulation yields a 
steady state value very close to this (the difference is due to the mechanical resistance load). Figure 17 
expands the Figure 16 time scale to show the three stator currents with 120° phase difference, daring 
motor start-up. 



Figure 18-20 shows the currents in the five rotor bars and the rotor velocity. Recall there exists 2n/5 
phase difference between currents in neighboring bars. This is clearly shown in Figure 19, which 
represents the motor starting moment. While the 60 Hz frequency of the stator currents generate a 
20 constant rotational velocity of the rotating magnetic field, the frequency of currents in the rotor bars 
decrease continuously as the rotor velocity increases. This is related to 'slip' in induction motors, the 
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nonnalized difference between the electrical angular velocity of the aii gap MMF established by the 
stator currents, and the electrical angular velocity of the rotor (Krause and Wasynczuk, 1989). Slip is 



defined as 



s = 



<D-(Q r 



— —1 



(22) 



CO. 



5 where an is the synchronous speed, or the speed of the stator currents, and ©r is the speed of the rotor. 
The magnitude and frequency of the currents and voltages of the rotor depend on the relative velocity 
between the rotating magnetic field and the rotor. In these simulations, this relative velocity rmxirnizes 
at t - 0, where the slip is unity. As the rotor velocity increases, the relative velocity and the slip 
decrease, suggesting that the decrease of amplitude and frequency of rotor bar currents in Figure 1 8 arc 

10 probably due to the decrease of slip. If o s - Or, slip s = 0 and no current is induced in me rotor bars 
(hence no torque). However, the. steady state currents of the rotor bars in Figure 18 are not zero, (even 
though there is no external load) because of the fractional load of the bearing modeled as a resistance 
R;c in Figure 1 5. If an external load is applied to the motor axis, the slip should increase and therefore 
the current and voltage in the rotor bars should also increase. Figure 20 shows the currents in the rotor 

15 bars during steady state. 

All simulation results shown above arc for a healthy motor. When rotor bars break, currents, 
velocity, and torque will deviate. Because we have a one-to-one correspondence between bond graph 
elements and machine components, it is possible to represent broken rotor bars by increasing the rotor 
bar resistance R,. In modern squirrel cage induction motors, bars and end-rings contact the rotor core. 

20 Due to this available current shunt, currents in a broken bar are not zero (Manolas and Tegopolous, 
1997), i.e., the resistance is not infinity. Figure 21 shows the stator currents and rotor velocity for a 
rotor with the third rotor bar broken. During the transient rise tune, the rotor velocity increases, and 
exhibits oscillations. Even at steady state, there exists periodic deviations of rotor velocity. With these 
deviations, the amplitude of the currents in the stator coils also change. These changes are more clearly 

25 presented in Figure 23; for comparison, a corresponding healthy machine simulation is shown in Figure 
22. Figure 24 plots the currents in each rotor bar, with bar 3 assumed broken. From Figure 24, the 
induced currents are largest in the two rotor bars nearest the broken bar. Figure 25 compares the torque 
characteristics of the healthy machine and broken bar machine. The rotor torque oscillates in the broken 
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bar machine, even at steady state. During startup, the osculation of torque is larger in the broken bar 

machine than the healthy machine. 

Simulations of an induction motor with a short circuited stator coil are shown m Figure Zo~Z8. 

In these simulations, the resistance of the snorted coil decreases, and me coil current, the magnetic 

5 fields, and the induced currents in the rotor bars also change. Figure 26 shows a difference in rise time 
of rotor velocity between the healthy machine and the stator coil short-crremted machine. Figure 27 
shows the rotor torque for both healthy and shorted rnachines. The overall trend of the torques are 
similar, but there exists small amplitude and relatively high frequency osculations in the short-circuit 
case. These oscillations are also seen in the rotor bar currents, Figure 28, compared with the rotor bar 

1 0 currents of the healthy machine, shown in Figure 16. 

A bond graph model of a squirrel cage induction motor was constructed, based on a prior bond 
graph by Ghosh and Bhadra (1993), that exhibited a one-to-one correspondence between the bond 
graph elements and real system components. Included were stator coil windings for three phases, 
mathematical transforations to incorporate two reaction theory, magnetic state variables to represent 

15 magnetic interactions between stator and rotor, individual rotor bars and attributions to the total rotor 
torque and velocity, and mechanical inertias and resistances. The srmulauons in this article had five 
rotor bars. Using this model, simulations of a healthy rmchine were compared to simulations of 
machines with a broken rotor bar breakage and a shorted stator coil. The degraded machine simulations 
predicted oscillations in currents and angular velocities, seen in real motors. 

20 Most induction motor designs employ three phase excitation of the stator. For a rotor with 

more bars, the bond graph of Figure 15 can be easily altered. More rotor bars can be included in Figure 
15 by adding additional pairs of power pathways to the right of the 2-port C's, such that n power 
pathways fan out from both a and p rotor phases. For the new value of «, these power pathways must 
update equations (11) and (12) governing moduli rm\ for the modulated transformers MTF:ro7* and 

25 equation (18) governing modulus i\ of the modulated gyrators MGY:r*. To update the 
electromechanical torque, in equation (17) we must replace the 5 in the upper index of the sum and the 
square root argument in the denominator with the new value of n. 

A Srr." nH Rvwmplary Appliffltinn 
30 This subsection refers to equations 1-6 in the detailed description. In addition, the naming 

equation numbers refer to equations vvithin this subsection. Further, an appendix is attached which is 
•referenced in this subsection. 

In a former embodiment of an induction motor, the bond graph model of a squirrel cage 
induction motor from above is adjusted. This model includes stator windings for 3 phases, two-reaction 
35 theory, magnetic interactions between stator and rotor, individual rotor bar contributions to rotor torque 
and velocity, mechanical inertias, and resistances and losses. Although this model does not include 
certain critical phenomena of the induction motor - e.g., magnetic field with rotor eccentricity or rotor 
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dynamics - this model is simple and can illustrate how to apply Shannon's communication theory to 
machine systems. 

In the system shown in Figure 29, MSe:V a , MSe:V b and MSe:V c indicate the 3-phase 
alternating voltage applied to the motor. The resistor element RJU models resistive losses in the stator 
5 windings of the motor. The gyrator GY:n 5 models the transition from the electric to the magnetic 
domain of the power flow in the system. The modulus of the gyrator n s equals the number of turns of the 
stator coil. The battery of transformers TF:mk convert the 3-phase into a rotating phasor vector. The 
two-port capacitance elements C represent the interaction between stator and rotor fields. 

In the rotor, electric voltage is induced in the metal bars by time varying flux cutting the bar 

1 0 circuits. This represented as the battery of gyrators, which have moduli n, related to the number of turns 
of the rotor. The modulated transformers MTTanrt relate angular position of the rotor relative to the 
flux field. The resistor elements R:R, represent resistive losses in the rotor circuits. Modulated gyrators 
MGYrric convert bar currents on the rotor bars into torque; this is the magneto-mechanical interaction. 
The moduli for these gyrators depend on fluxes in the rotor. The final transformer TF:mm is related to 

15 the number of magnetic poles in the system. Power lost by bearing friction is accounted for by the 
resistance R:c. The remainder of the power drives the output shaft 
Moduli in the bond graph (Figure 29) are given as follows [6]: 
1) Moduli of three phases are (jnt) : 

iii 

20 2) Constitutive laws for two port C fields arc : 
where, * = CL, fi and 



L,L,-L m " " L s L r -Lm 2 ' ' L s L r -L* ' 
L s is stator self inductance, L m is mutual inductance and^, is rotor self inductance. The gyrator 
25 moduli is the number of stator coil turns, and gyrator moduli n r is the number of rotor coil turns. 
3) The modulated transformers MTFrmrk are : 

&+ 2(k-\,„ , k = lt 2,..., n 



'.I— cos< t> + — > 

Vn [ 7i J 



12 
in 

V« [ n J 



(9) 



where n is the total number of bars. 
30 4) The moduli of the modulated gyrators MGY:rjc are : 
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5) The modulus for transformer TFrn^, is : 
P p 

m m = — : P p is number of poles 



(ID 



State equations were derived from the bond graph (figure 29) with n s j = - n s2 - n rI - n r tcr- n r . 
The results of state equations are 



. 1 h 
5 0 o =— 7 



(12) 



Par = -{ mr S R r\ +»W 1 1 R r2 +mr 3 2 ^ r3 + wr 4 2j ^4 +7Kr 5"^5 

_ A- 

-(mrjinr,,^,., + mr,mr,R, 2 + mr 3 mr t R r , +mr A mr $ R r4 + mr s mr w R r} >— 7- 

, 1 h 

-(mr, r, +mr 2 -r 2 +mr 3 r 3 +mr 4 r A +mr s r s ) 



(13) 



1 



10 



- (mr 6 • r x + mr 7 • r 2 + mr z • r 3 + mr 9 • r 4 + mr l0 • r 5 )^ ^ — 



(14) 



1 — 



n s l^fl 1 jtf 



2 fl 2 



2 r *j 



2 \ 



1 2 n 2 

V 72 1 J 



1- ^ 



4 J) 



(15) 



1 — 



R. 



2 > 



1 — ^ ^ 

5' n *(Jv n >~ +R * ) 



' n s 2 +R sl \m 3 



m, 



s J) 



(16) . 



/i = (r, mr, +r 2 mr 2 + r 3 • mr } + r, • mr A +r s mr s ) — — M„. 

+(r, -mr 6 + r, • mr 7 + r 3 • mr 8 + r 4 • *nr 9 + r 5 ■ mr 10 J M A — -c 



(17) 
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where the magnetic state variables are rotor angular position 0 O and momentum h and stator and rotor 

phase fluxes , , <P*r 9 fir • 

Simulation of a squirrel cage induction motor employed MATLAB®'s Runge-Kutta 4 th order 
method with a time step 8t = 10~ 5 seconds. Values of the system parameters presented in Table 2 

5 were given by Kim and Bryant [6,7]. 

Using this model simulations were performed for an ideal machine, which has no faults and 
functions perfectly according to designer's specifications, and a degraded inachine. The ideal machine 
will serve as a reference of desired dynamic behavior. The degraded motor will exhibit common 
degradation modes, including rotor bar breakage and stator coil snorts. We will excite the ideal and 
1 0 degraded machine models with identical test signals, record these signals, and then estimate the noise as 
the difference between degraded and ideal machine responses to the same test signal. 



Table 2 System parameters of a two pole, three phase squirrel cage induction motor. 



Parameter 


Description 


Ideal 


Degraded 


Rsa> Rib, Rk 


Stator coil resistance [CI] 


0.0788 


0.00079 - 0.0709 


R,J t ... JtrtO 


Rotor bar resistance [CI] 


0.0408 


0.0412-4.0800 




Input voltage arnplitude [V] 


230 




f 


input voltage frequency [Hz] 


60 






Stator Reluctance [l/H] 


5.85x10* 




% 


Rotor Reluctance [l/H] 


563 






Mutual Reluctance [l/H] 


5.41x10* 




n 


Number of rotor bars 


5 




n s 


Number of effective stator coil turns 


100 


1-90 


n r 


Number of effective rotor coil turns (bar) 


1 




c 


Mechanical resistance [N-s/m] 


0.15 




J 


Mechanical inertia [kg.m 2 ] 


0.1 




m n (-P/2) 


Number of poles □ 2 


2 





15 

Figure 30-32 show sample simulation results for a nominal or ideal motor, Le., a motor 
without faults. These simulations arose from the model of equations (7) to (17), with the "Ideal" 
parameter values of Table 2. Plotted are selected motor state variables versus time, beginning with 
motor startup, i.e. the moior voltages were switched "on". The rotor velocity rises slowly to a steady 
20 state value of about 377 rad/sec, as the momentum and all other state variables reach steady state. 
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Theoretically, when 6OH2 alternatuig inputs are given to a two pole AC motor, the output velocity 
should he close to 3600RPM (= 377 rad/sec). Rotor and stator fluxes of a phases are shown in figure 
3 1 and 32; the P fluxes are similar. Flux amplitudes increase to steady state, consistent with the 
angular velocity. 

5 Various feults can he developed in motors. For example, stator coil shorts cause overheating, 

increasing core losses [8]; rotor bar breaks or cracks in the die-cast rotors cause very large electrical 
resistance [6, 7, 9]; and bent or cracked shafts make the rotation wobble [10]. 

In this article, we will focus on a broken rotor bar, and shorted stator coils. When rotor bars 
break, steady state velocity and torque of the rotor will deviate from the ideal response. With the bond 

10 graph shown in Figure 29, a broken bar can be incorporated into the model by increasing selected rotor 
bar resistances The range of deviation is given in Table 2, last column. Figure 33 (upper curve) 
shows the response of the motor with a broken bar after being switched W. In this case, the rotor bar 
resistance was increased 10 times from its nominal value of 0.0408 0,ohrns, toR,= 0.408 D.ohrns. 
Plotted is the angular velocity versus time, from start up. Here the angular velocity increases, during a 

15 transient time characterized by deviations of rotor velocity. 

Figure 33 (bottom curve) shows the simulated startup (step) response for a motor with a rotor 
bar having resistance increased 100 times, to R, = 4.08 D-ohms. This curve shows increased and 
persistent oscillations, compared to Figure 30 for the ideal machine. 

It is well established that when rotor faults occur, rotor harmonic fluxes are produced which 

20 induce currents in the stator at frequencies of f[k/(P p /2).(l-5>± s\. Here /is the supply 
frequency, P p is the number of poles , k =1,2,3,4. . and s is slip defined as [11, 12] 
S ^^LZ^L (18) 

In equation (18) ® 3 is the synchronous speed derived from the frequency of the stator currents, and 

<x> r is the angular speed of the rotor. Slip can have a value from 0 to 1 . 
25 Figure 34 (a) ; constructed by applying a Fourier transform to the steady state portion of the 

simulation results of Figure 33 (bottom curve), shows some frequencies of the stator current of phase A 

in the vicinity of the excitation frequency (60Hz). Figure 34 (b), scanned from reference [13], shows 

spectral densities of typical currents versus frequency measured from a motor with three broken bars. 

Comparison of Figure 34 (a) and (b) show similar shape and location of spectral peaks. 
30 Figure 35 compares the torque characteristic of the ideal machine and broken bar machine R r = 

0.408 □.ohms. The rotor torque oscillates whenever the rotor velocity oscillates. Due to the rotor 

asymmetry the level of pulsating torque is increased [6, 12]. 

The average power in a signal *(0 , of duration T can be estimated as [14] 

1 0 
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or as 



If x n is a sequence sanded from x(t) at equally spaced discrete instants. The power spectral density, 
the magnitude squared of the Fourier transform of *(0 , is given by 
5 S(f)=\X(ff ™ 

where 

—aj 

For discrete x„, we employed a fast Fourier transform to obtain Ai. 

The total power can be calculated in the frequency domain, or in the rime domain by Parseval's 
theorem [14] 

I0 \x{tfdt=)\X(jfdf (22) 
The discrete form of ParsevaTs theorem is defined as [14] 

iw 2 4ikr- . (23) 

Equation (5) can be rewritten as 
15 C = )§Og 2 (S + lf>-)og 2 N}d0. (24) 

0 

Combination of the original signal power spectral density (5) and the noise power spectral density (AO 
represents the signal power spectral density 

S.=S + N < 25 > 

from the uegraded machine. 

20 Shannon [1] assumed a Gaussian white noise statistically independent of the signal. To remove 

this restriction, we need to calculate the noise power directly from the time domain signals. In the time 
domain, the noise is defined as the difference between actual and ideal signals 
n(t) =x(t)-x 0 (t). ( 26 > 
Here x(t) is the output of the degraded machine, and x o {i) is the ourput of the ideal machine. As 

25 demonstrated in the Appendix for this subsection, removal of the independence restrictions between 
signal and noise admits negative channel capacities. 

Power spectral densities £ and N' can be defined as the magnitude squared of the Fourier 
transforms for signal x(t) and noise "(0 respectively. To calculate the channel capacity with these 
values, we must replace {S + N) in equation (24) with SI , and N with N' , to have 
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Figure 36 (upper line) shows the power spectrum of the rotor velocity of the ideal machine as 
shown in Figure 30 (npper line), and defined in section 4.1 . This figure was constructed by applying a 
fasi Fourier transform to the angular velocity data of Figure 30 (upper line). In this case, we assume 

5 zero noise, and thus the system functions perfectly, according to the design specifications. Using 
equation (27), we get an infinite channel capacity for an ideal machine system, since by definition, me 
noise and noise power are zero. 

If there are faults such as broken rotor bars as mentioned in section 4.2earlier, the power 
spectrum will change as noise contaminates the signal. Figure 37 shows the startup response x(t) of 

10 the machine with a broken bar (upper line), and the noise in the time domain from the degraded 
machine, defined by equation (26). This noise (lower curve and magnified in Figure 37) is the 
difference between the degraded nnchine's response curve in Figure 37, and the startup response of me 
ideal machine in Figure 30. The presence of several frequencies is evident. Figure 36 (dots in the upper 
line) shows the power spectra (signals and noise) of a degraded machine, with a cracked rotor bar, Note 

15 that the power spectrum of the degraded machine signal x(t) nearly overlaps the power spectrum of 
the ideal machine signal x c (t) ; in the figure, the two almost coincide. In our model, we increased the 
resistance of broken (cracked) bar by 10%. Using equation (27), we obtained a channel capacity of 
1.1 xlO 6 (bits per second). Here the integration bandwidth CO m equation (27) was equated to the 
entire sampling bandwidth (5xl0 4 ifc) based on the Nyquist's sampling rate, where St = 10 

20 seconds was me time step employed in the numerical solution routine, m this procedure, we viewed the 
numerical solution's data points as a "sampled" signal, with sampling interval equal to the numerical 
method's time step. The Nyquist's sampling rate gives the smallest bandwidth associated with the 
sampling interval St = \0~ 5 seconds. 

In equation (6) for entropy rate R, S i represents the average power of the output signal from 

25 the healthy ruachine and N t represents the largest acceptable deviation, i.e., a tolerance on the noise. 
The signal bandwidth (#,-) was equated to G>\ see the previous paragraph for justification. The 
functional requirements of the machine tetennine the noise or error tolerance N { demanded by the 
inacbine to work satisfectorily. For example, if we have an application wherein the rmximum allowed 
error or tolerance must be within 10% of the signal of the ideal machine, and if we employ the same 

30 bandwidth as forme channel capacity, then from the equations (6) and(19), the information rate (R) is 



7? = <» i log 2 (S i /#,.) = 5xl0 4 log : 



U0.\xx(tfdt 



Vt 



= 3.3xl0 5 (bits per second) (28) 
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With a channel capacity of 1 • 1 * 1 0 6 (bits per second), the result for R above satisfies the condition of 
R^C . This suggests a still "healthy" machine. 

If the resistance of the broken bar increases 21 times to R, = 0.8568D, then the channel 
capacity drops to 2.4xl0 5 (bits per second) below the required (R) of 3.3x1 0 5 (bits per second). 
5 Since this result doesn't satisfy R^C, according to Shannon's theorem, the machine is 
malfuncuonrog. As the magnitude of the fault (bar resistance) increases, the channel capacity 
diminishes. 

Figures 38(a)-38(d), show selected power spectral densities of the stator current of phase A at 
steady state, for selected bar resistances. These figures are similar to Figure 34 (a) and are often used as 

10 diagnostic indicators. Side bands are absent for smaller values of but start to appear after tbe rotor 
bar resistance equals approximately 0.7670Dohms (1780%). Figure 38(e) also plots the channel 
capacities versus the percent change of bar resistance from the ideal value given by table 2. The dashed 
line indicates the 10% noise power tolerance {R = 3.3 x 1 0 5 ) estimated in the previous paragraph. Here 
percent change is defined as %R bar = (R bar -R 0 ^ , where R^ is the current value, and 

15 Rfr r the ideal. The channel capacity of point (d) in figure 38(e) has a negative value; see the Appendix 
for reasons. As shown in figures 38(c) and 38(d), significant side bands with large intensity begin to 
appear in the power spectra, wherever the channel capacity curve sinks below the 10% information rate 
line (dashed curve). From a practical standpoint, for industrial grade rugged machines such as motors, 
we would begin to notice errors when these exceed 10% or more in the motor's output velocity. Thus, 

20 1 0% was chosen as the critical velocity condition. 

The curve of figure 38(e) can be separated into regions with three distinct slopes: region 0, 
which connects the infinite channel capacity of the ideal system to that of "real" systems; region I, with 
stable C and "healthy" operation (region I would be associated with the normal life cycle operating 
region of the system); and region U, where C declines to the (dashed) failure line. Note that the marked 

25 change in the slope of C or the rapidly dinnmshmg values of C, going from region I to ^ could presage 
failure. Figures 38(b)-(d) suggest that once side bands appear, the slope of C becomes noticeably more 
negative. 

Signal based diagnostic methods, that trigger upon detection of side bands, at earliest would 
notice the broken bar fault at point (b) in figure 38(e); figure 38(b) suggests that detection of the tiny 
30 side band would be difficult. In contrast, the channel capacity curve's knee - where the slope abruptly 
changes - occurs at 1500%, before 1780% of figure 38(b). Here the abrupt change in slope might be 
easier to detect 

Simulation of an induction motor with short circuited turns on its stator coil is shown in figure 
39. Here stator resistance R«i of phase A was decreased 50%, and the effective number of turns 
35 represented by gyrator modulus n s was similarly decreased from 100 to 50. In this simulation, only one 
of the stator coils has shorted turns. In the model and physically, as turns are short circuited, the 
resistance in that coil decreases, and the effective number of turns also decreases. Figure 40 shows the 



013270.00015:124601.01 



WO 02/27418 



PCT/US01/29935 



35 

various signal and noise power spectra of this shorted machine. Note again that the power spectra of the 
ideal and degraded machines nearly coincide. The shorted coil seems to affect the angular velocity 
relatively less than a broken bar. The two startup responses in the figure 39 are nearly the same, 
especially at steady state. Essential differences at steady state are in the signal's phase, generally not 
5 contained in power spectra. By taking the difference between degraded and ideal response, the 
information on phase differences is conveyed in the noise power spectrum, in addition to die magnitude 
information. 

When stator coil turns short out, we observe only a rise in some of the frequency components 
which already exist in the stator current spectra of an ideal machine [15]. Figures 41 (a) and (b) shows 
1 0 spectral content of the steady state stator current of phase A, from simulauons of the bond graph model 
For comparison, spectra from Gojko and Penman's model [15] are also shown as figures 41 (c) and (d). 
Figure 42 shows spectral content of stator currents for two shorted coils, phases A and B. In figures 43 
and 44 are plotted the channel capacities versus percent change in the coil resistance, for shorting of 
phases A, and A and B, respectively. Again information rate for the 10% noise level on angular velocity 

1 5 is shown as the dashed line. 

Similar to the broken bar case of figure 38, figure 43 exhibits a ' Wfhy" region I, with stable 
channel capacity, and a region II with sharply dimrmshing channel capacity. Again the sharply changed 
slope of region II could prognose failure. 

In this article, it was demonstrated how Shannon's theory of connmmication could be applied 

20 to machinery, to utilize Shannon's powerful theorems. Concepts of rate of information rate and channel 
capacity were reviewed, and applied to an mcucnon motor. At the heart of the method is machine 
'noise", estimated as the difference between actual and ideal responses. By subtracting the ideal 
■ response from the actual response *(0 , the noise signal contains only information about the 

faults. From the noise and signal were calculated power spectra, used in equation (27) for channel 

25 capacity. Rate R from equation (28) serves as a critical values dependent on the system's tolerance to 
errors, here called "noise". The channel capacity was calculated for a motor with shorted stator coils 
and broken bars. The channel concept agreed with other existing fault monitoring methods, but results 
suggest that it could detect faults much earlier. It can be concluded from this study that the channel 
capacity concept could serve as an effective tacriinrnator of motor and machine system health. 

30 It is to be understood that the forms of the invention shown and described herein are to be 

taken as the presently preferred embodiments. Elements and materials may be substituted for those 
illustrated and described herein, parts and processes may be reversed, and certain features of the 
invention may be utilized independently, all as would be apparent to one, skilled in the art after having 
the benefit of this description of the invention. Changes may be made in the elements described herein 

35 without departing from the spirit and scope of the invention as described in the following claims. 

As such, a method for diagnosing the state of a system is described, in view of the 
above detailed description of the present invention and associated drawings, other 
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modifications and variations will now become apparent to those skilled in trie art. It should 
also be apparent that such other modifications and variations may be effected without 
departing from the spirit and scope of the present invention as set forth in the claims which 
follow. 
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Nomenclature 

C channel capacity 

c viscous resistance coefficient 

f supply frequency 

5 H information entropy 

h angular momentum [N-m-sec] 

/ moment of inertia 

Ls, Lm, It stator self inductance, mutual inductance and rotor self inductance 

M , M p a and |3 axis magneto motive force 

10 mi~ms moduli of three phases 

ji modulus of gyrator (number of coil turns) 

n s number of effective stator coil turn 

ttr number of effective rotor coil turn 

n(t ) noise in time domain 

IS P p number of poles 

p probability of occurrence 

R entropy rate 

jfc, Rr stator and rotor resistances 

r ,~r 5 modulated gyrator moduli of rotor 

20 SK,, ^m, stator reluctance, irnimal rehiaance, rotor reluctance 

S 9 N average power of the signal and noise 

S„ average power of the signal including noises 

s slip 

Va. Vb, V c sinusoidal input voltages 
25 x(t) output of the degraded m."<*mc in time domain 

x o (t) output of the ideally healthy machine in time domain 

CO bandwidth 

P e > Pp a and p axis fluxes 

□ flux linkage 

30 Vo,, V ps a and p axis stator voltages 

i^, i^j a and p axis stator currents 

i^, ip r a and p axis rotor currents 

Rs, Rj stator and rotor resistances 

Lj, Lm, Lr stator self inductance, mutual inductance and rotor self inductance 

35 T c . Tl electro-magnetic torque and mechanical load torque 

J moment of inertia 
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c viscous resistance coefficient 

^ ^ electrical and mechanical angular velocities of the rotoi 

P number of pole pairs 

X flux linkage 

5 mi-m 5 moduli of transformers for 3 phase to 2 phase transformation 

V a ,V b ,V b sinusoidal input voltages 

id-ms transformer moduli 

1^ current in the k A rotor bar, 

m magnitude modulus that depends on the total number of bars 

I o n modulus of gyrator (number of coil turns) 

$ magnetic flux [Weber (Wb)] 

M magneto motive force [Ampere (A)] 

p permeance of circuit element [Henry (H)] 

R reluctance of circuit element [1/Henry (H* 1 )] 

15 ei,e 2 effort 

fi,f 2 flow 

n 3 number of effective stator coil turn 

ni number of effective rotor coil turn 

h : artg^ ar momentum [N*m-s= : kg"m 2 -sec] 

20 Rs_arpha, 
Rs_beta, 
R«,Rsa, 

Rsb, Rsc, R electrical resistances 



25 
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APPENDIX FOR FIRST SQUIRREL CAGE INDUCTION MOTOR EXAMPLE 

In this section, for a rotor with n bars : we prove 



(A.1) 

5 From equations (15) and (16), the transformation matrix tones its transpose is 

4/0 



A T A = m 



sin0 sin 
cos# 

cos 



cos^ cos^+— j cos^+— j ... cos[0 + ^- 

(a 2 A 
: 6 + — 

I n ) 



sinj 0 + — 



1 2(n-l)^ j 



(A.2) 

Hie result of the multiplication is a square matrix of dimension 2, 



A T A = m 2 



^21 ^22, 



10 (A3) 
where 



^ -g^»{*+at=a=}. cm» 

Equations (A.4) can be rewritten using double angle trigonometric formulas: 



15 ill-Z 



n i + cos {2 e^} s in -i 



= ?-tZcos{20+jf } 



k=l 



"2 '2. 



j=0 



~k=l j=o 



822= I" 



l-cos{29+ v \ } 



n-1 



4tc 

- 2 - 2 2cos{2e+j-) 



(A.5) 



Via formulas 1.341-1 and 1.341-3 in Gradshtcyn and Ryzhik (1980), the sums of sine and cosine terms 
on the right sides of equations (A.5) are zero, for n > 3 . Thus 
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S„ = 



1 + icos2/9 

2 2 

1 + cos 20 

w 

2 

(A.6) 

i-I C os20 
2 2 
1- cos 20 
n 

2 



;n = l 
;n = 2 
;n>3 



;n = l 
;n = 2 
;n>3 



(A.7) 



5 15 — S 7 . — 



-sin 20 
2 

sin 26 



;n = l 
;n = 2 
;nS3 



(A.8) 

Therefore, we can conclude 



A T A = m : 



?2l S 12. 



- m' 



i 0 
2 

0 £ 
2. 



(A.9) 

10 for the rotor with more than 2 rotor bars, i.e., n > 3 . 
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Appendix For Second Example of Induction Motor 

Range of Channel Capacity Values 

In equation (27), average power of the output signal including noise is defined as 



c 



(A.1) 

Using equations (3) and (4), we get 
2 r 

s: = S+ AT' +-Jk(0 •»(')]*. 

■* 0 

(A.2) 

10 Shannon assumed statistical independence of X B (t) and »(0 , which made the last term of equation 
(A.2) vanish. In this article, we will remove this restriction, allowing for forms of noise 
n (/) = -K.v,(/).(0<K<l) 

(A.3) 

that can extinguish the signal, such that 

15 x(O = *.(0+«(') = (l-Kk(0 
(A.4) 

Here the amplitude (l " K) of x(t) diminishes with increasing noise: this can affect the signal power 
significantly. With the values just derived, 

. C-xyijkotf* f ,y 
s - 1 0 -= i-- 



N' „ 3 l r rr V 



20 (A.5) 

can be less than unity, and the channel capacity can have negative values: as noise power proportional 
to K- 2 increases, the output power proportional to (l-K) decreases. 
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We claim: 

1 . A method of modeling a mechanical system comprising a plurality of physical components, 
comprising: 

preparing a model of a mechanical system in which at least a portion of the physical 
5 components of the mechanical system are individually modeled, wherein the model is 

configured to output data representative of the condition of the mechanical system in response 
to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predetermined operating 
conditions during use; 

1 0 modifying the model such that the outputtsd data of the model in response to the 

predeteimined operating conditions is representative of the condition of the mechanical system 
in response to the predeteimined operating conditions. 

2. The method of claim 1; further cornprising predicting a failure time of the mechanical system 
using the modified modeL 

15 3. The method of claim 1, wherein all the physical components of the mechanical system are 
individually modeled. 

4. The method of claim 1, wherein a possible fault for each of the mdividually modeled physical 
components are incorporated into the model 

5. Thcmethod of claim 1, wherein a potential failure for each of the individually modeled 
20 physical components are incorporated into the model 

6. The method of claim 1 , wherein a plurality of possible failures for the individually modeled 
physical components may interact, rendering failures not specifically associated with any single 
component, but arising from interactions between components. 

7. The method of claim 1, wherein the condition of the modeled mechanical system is 
25 represented within the model as noise. 

S . The method of claim 1 , wherein the condition of the modeled mechanical system is 
represented within the model as noise, and wherein the condition of the modeled mechanical system is 
determined by calculating a signal to noise ratio for the model. 

9. The method of claim 1 , further comprising: 
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calculating the channel capacity of the modeled mechanical system, wherein the channel 
capacity is representative of the design of the system and the present condition of the 
mechanical system; 

calculating a rate of information for a predetermined job to be performedby the mechanical 
5 system, wherein the rate of information is representative of the speed, loads, complexity and 

desired accuracy of the job; and 

comparing the rate of information to the channel capacity, wherein if the rate of information is 
less than or equal to the channel capacity the model will output data indicating that the 
mechanical system is capable of performing the job at the appropriate speed, load and 
10 accuracy. 

10. A computer implemented method of modeling a mechanical system comprising a plurality of 
physical components, the method comprising: 

preparing a model of a mechanical system in which at least a portion of the physical 
components of the mechanical system are individually modeled, wherein the model is 
15 configured to output data representative of the condition of the mechanical system in response 

to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to pretetermined operating 
conditions during use; 

modifying the model such mat the outputted data of the model in response to the 
20 predetennined operating conditions is representative of the condition of the mechanical system 

in response to me predetermined operating conditions. 

11. A carrier medium comprising computer instructions, wherein the program instructions are 
computer-executable to implement a method of modeling a mechanical system comprising a plurality of 
physical components, the method comprising: 

25 preparing a model of a mechanical system in which at least a portion of the physical 

components of the mechanical system are individually modeled, wherein the model is 
configured to output data representative of the condition of the mechanical system in response 
to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predetermined operating 
30 conditions during use; 
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modifying the model such that the outputted data of the model in response to the 
predetermined operating conditions is representative of the condition of the mechanical system 
in response to the predetermined operating conditions. 



5 12. A method for diagnosing a state of a system, the method comprising: 
measuring a signal from the system; 

comparing the signal from the system and an expected signal to determine a noise signal 
associated with the signal from the system; 

determining a signal strength associated with the signal from the system; 

1 0 determining a rate of information, the rate of information associated with a desired operability 

of the system; 

detennming a channel capacity from the noise signal and the signal strength, the channel 
capacity being a function of a frequency spectrum of the signal; 

comparing the rate of rnformation to the channel capacity to diagnosis the state of the system 

15 

13. The method of Claim 12 wherein the expected signal is a signal measured from an exemplary 
system operating in a known state. 

14. The method of Claim 12 wherein the expected signal is the output of amodel. 

20 

15. The method of Claim 14 wherein the output of the model is adapted to approximate the measured 
signal. 



16. The method of Claim 12, the method further comprising: 
25 repeating the steps of the method over time to determine a set of diagnoses, 

detennming a prognosis of the system from the set of diagnoses. 
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